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APPLICATION OF DESCRIBING- FUNCTION ANALYSIS TO THE 
STUDY OF AN ON-OFF REACTION- CONTROL SYSTEM 
By Edgar C. Lineberry, Jr. , and Edwin C. Foudriat 

SUMMARY 


An analytical study was made of an automatic reaction- control sys- 
tem for the upper stages of a missile to determine limit cycle charac- 
teristics, corresponding duty cycles, and the effects of the various 
system parameters on these quantities. A nonlinear servo analysis 
(describing function) technique was used to obtain a mathematical repre- 
sentation of the nonlinear components of the system. 

The results obtained by this analysis are compared with the results 
obtained from an analog simulation including the jet reaction- control 
hardware. The good agreement between the results of the two studies 
tends to indicate the feasibility of using such an analysis technique 
in the early design phase of a reaction- control system, since by so 
doing, the parameters for good system performance can be determined 
quite readily. 


INTRODUCTION 


Reaction- control- system characteristics and requirements are usually 
determined by simulating the system and the dynamics of the vehicle on 
an analog computer. The system parameters for final design are obtained 
from the tests. Such techniques, while necessary, are undesirably time 
consuming especially for the early design stages when quick results and 
a clear-cut picture of the effects of various system parameters are 
required. Therefore, there is a need for an analytical technique which 
demonstrates the effects of the various system parameters and gives a 
good approximation to the actual results and the system modes of operation. 

The analysis of a closed- loop system in which all of the system 
components can be represented by linear differential equations can be 
made by using conventional methods. However, if nonlinear components 
(components which display discontinuities in their operation) exist in 
the system, the techniques must be revised in order to handle the non- 
linearities. Such is the case with the react ion- control system studied 
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in this paper. A theoretical analysis is made of the reaction- control 
system to determine the possible stable and unstable modes of opera- 
tion, their corresponding duty cycles, and the effects of the various 
system parameters on these quantities. The describing- function tech- 
nique (refs. 1 and 2) is used to obtain a mathematical representation 
of the nonlinear components of the system. 

It is the purpose of this paper to indicate the validity of an 
analysis with the use of the describing- function technique and the 
feasibility of using such an analysis in the early design phase of a 
reaction- control system. Inherent in the technique is the capability 
of presenting clearly the effects of the system parameters upon the 
system operation. 


SYMBOLS 


bm 

D 

F j 

% 

G 

I 

J 

k 

K e 

K e 

K a 


Fourier series coefficients 
reaction- control- thrust moment arm, ft 
misalinement- thrust moment arm, ft 
control- system dead zone, deg 
reaction- control thrust, lb 

misalinement thrust, lb 
linear gain 

moment of inertia, slug-ft 2 
vector operator, vrr 

ratio of valve turn-off to turn-on time 

position gain, deg/deg 

. deg/deg 
rate gain, — — — 

deg/sec^ 

acceleration constant, 


deg 
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N nonlinear gain 

R ratio of misalinement moment to control moment 

s Laplacian operator 

t time, sec 

X switching system input 

Y switching system output 

a, 0 nonlinear switching angles, radians 

£ damping ratio 

9 pitch angle, radians 

T) angular difference between on and off times, radians 

t valve response time constant, sec 

0 phase angle, deg 

cd frequency, radians/sec 

CDn undamped natural frequency, radians/sec 

CDf limit-cycle frequency, radians/sec 

Subscripts: 

1 control- system input 

j reaction jet 

r control- system output 

1 fundamental component 

ss steady state 

A dot over a symbol designates the time derivative of the variable. 
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SYSTEM TO BE ANALYZED 


The system being investigated controls the attitude of a missile 
by producing moments to overcome disturbing moments encountered by the 
vehicle. The restoring moments are produced by small reaction jets 
mounted near the base of the vehicle with thrust vectors normal to the 
vehicle longitudinal body axis. In this analysis the system is assumed 
to control the pitch attitude of the upper stages of a missile in the 
portion of its controlled flight where aerodynamic forces and moments 
can be neglected. The remaining moments consist of thrust misalinement 
resulting from the fixed main engine and the reaction- control moments. 
With the assumption that the vehicle is rigid, a single-degree-of- 
freedom system can be assumed, and the equation of motion becomes 

10 r = F^bj + -^m^m (^- a ) 

and after burnout of the main engine is 

I0 r = Fjbj (lb) 

A block diagram of the control system is shown in figure l(a), 
showing the nonlinear components, a reaction motor, the vehicle dynamics, 
and the feedback system. The feedback system includes both perfect 
rate and position feedback and may be obtained physical! y from a rate 
and attitude gyro combination. 

The nonlinear components consist of an electronic relay with a 
threshold sensitivity or dead zone ±D and a two- stage electro- 
pneumatic valve system. The operation of the system can best be illus- 
trated by a time sequence of operation as shown in figure l(b). When 
the error signal exceeds the dead zone, as shown in the top trace, 
the electronic switch is actuated. A finite time t then elapses 
before the pneumatic- control valve is actuated and fuel flows out of 
the valve. Since the valve system turn-off time may not be the same 
as the turn-on time, the constant k, the ratio of turn-off to turn-on 
time, is used so that the turn-off time is kT. 


In the operation of the react ion- motor- control system, an addi- 
tional time lag occurs between the time the motor is turned on and the 
steady- state thrust is reached, due to the buildup of gas pressure in 
the chamber. This lag has been approximated by a second-order system 
whose transfer function is 


f j 


( F j)ss 




2£s 

7r + — — + 1 

2 % 


( 2 ) 
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A comparison between the approximated and the actual control force is 
shown at the bottom of figure l(b). 


DESCRIBING- FUNCTION TECHNIQUE 


For a theoretical analysis using the describing- function technique 
(refs. 1 and 2 ), the control system is divided into linear and nonlinear 
subsystems. Consider the closed loop system in the form shown in fig- 
ure 2 , with the block N representing the nonlinear components and 
G(s) representing the transfer function for the linear components. 

For the reaction- control system under consideration, the linear trans- 
fer function becomes: 


G(s) 


+ l) 




(3) 


In order for the describing- function analysis to accurately 
represent the operation of the system, it must be true that the linear 
sections of the system act as low pass filters; that is, the harmonics 
which may be generated in the nonlinear system are attenuated to a 
greater degree than the fundamental component of frequency. This tech- 
nique was applied in the range where equation (3) fulfills this 
requirement. 

For a limit cycle or steady- state oscillation to be maintained 
in the system, the frequency must be such that the system components 
produce a total phase shift of l80° when the gain in the system is 
unity. Using this relationship 

NG(Jod) = -1 (4) 

or 

G(jm) = - | (5) 

The describing function is used in the analysis to represent N. 

In order to determine N, an input of the form: 

X = Xj cos cut (6) 

is imposed at the input to the nonlinear block of the system. An out- 
put is assumed in the form of a Fourier series 
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Y = Aq + 2_, (A n cos nu>t + % sin nut) (7) 

n=l 

with the coefficients as functions of the nonlinear system N. Because 
of the assumption that the system acts as a low pass filter, the funda- 
mental component will be the only one considered to affect the system. 
Thus, the nonlinear gain N is the ratio of the fundamental component 
to the amplitude of the sinusoidal input, or 



N is deter mi ned by the describing- function method for the following 
cases: 

Case I. Describing Function with Reaction Control (Missile Coasting) 

In order to determine the gain of the nonlinear components composed 
of a dead zone ±D and a fuel-control valve with a time lag t, an 
input of the form 

X = X-l cos ojt 

is imposed on the system. From the previous discussion, the output is 

Yi = Aq + A^ cos cut + B-l sin tot 

with the coefficients Aq, Aq_, and B-^ to be determined. 

The system operation is shown graphically in figure 3 (a). The 
input triggers the electronic switch at -oq, the time t elapses 

(corresponding to an angle orr) before the valve opens at -a^, and 
fuel begins to flow. The turn off is at oq and 0 , 3 /, respectively. 

The identical operation takes place on the negative half cycle. Thus, 
for the positive half cycle, the angles at which the valve turn on and 
off are given by: 

-a-^ = -Oq + our (9) 

and 

a l' “ Oq + knur (lO) 

The phase shift for the nonlinear system is given by: 
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0t l + ( a i) CUT/., , N A 

02 = —(1 + k ) = 0 


It is possible to shift to a new axis, 02 shown in figure 3 (b), 

in order to simplify the evaluation of the coefficients of the Fourier 
series. The input now takes the form 


X = X^ cos (cut + 0 ) 


The output takes the form 


Y^_ = Aq + cos cut 


(12a) 


(12b) 


Since the output reaction thrusts are symmetrical and equal about the 
abscissa over the interval from 0 to 2jc, the term Aq is equal to 

zero. Bie A-^ term is given by 


A^ = ^ f f(o>t)cos cut d(cut) - ^ P f (cut) cos cut d(cut) (13) 
J 0 J $2 

where f(cut), the valve output, is represented by unity. Since 


f (cut) = 1 
f (cut) = -1 


(0 ^ cut te P^ 
^ cut % Jt) 


P 2 = * - 3 q 

substituting into equation (13) gives 


2 p ^1 2 P ^ 

^ / cos cut d(cut) - — / cos cut d(cut) 

^ 0 J it- Pi 


>1 - ai ' ~ g ( '° l) - «0 * f(k - 1) 
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and 


D = X]_ cos clq 


(18) 


or 


-1 D 

a 0 " COS Xi 


(19) 


Therefore, 


Pi = cos' 1 A + 2£<k - 1) 


( 20 ) 


Substituting equation (20) in equation (l6) and integrating gives 


A 1 ' ; sin^cos -1 £ + f(k - 1)] 


( 21 ) 


The describing function of the nonlinear components is given by 


N = 3^ = rtX^ Sln 008-1 + ^ 


Xl 2 


-00 


( 22 ) 


where 0 is given by equation (ll). 

The expression for N can now be substituted into equation (5)> 
giving 


±- sinfcos- 1 £♦£<*- l)]e J2 


-J<K<k+l) 


= G(j<c) 


(23) 


The values of <d and — which satisfy this relationship can be 

X 1 

determined graphically from the intersections of the curves G(jo>) and 
- ^ in a plot of gain against phase angle. Only one curve results for 
the linear expression G( jcc) for values of oo. For the nonlinear 
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expression, however, a different curve with — as a parameter, results 

X 1 

for each value of <x>. By observing, however, that the phase angle for 
the nonlinear expression N varies only with cd, a value of o> can be 
determined where the phase angles of the "two expressions G(jo>) and 

- ^ are equal. Thus, at this particular value of cn, the N and the 

G curves represent a possible intersection point. Rearranging equa- 
tion (23) and substituting the particular value of o> results in 


1 . r -1 d 

— sin cos — 

x i L x i 


f(k- 


1) 


jt 

4g 


Letting 


and 


C = = Constant 

4g 


t} = - l) = Constant 


(2L) 


(25) 


(26) 


gives 


r * 

— sin cos - ^ — + q = C 

X ! |_ X 1 


Applying trigonometric identities. 



tj + D sin 


= C 


(27) 


(28) 


Rearranging and solving for 


by letting 

X 1 


o 

B = 2CD sin tj + cos 


(29) 
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gives 


Thus two values of 
conditions. 


D_ 




1/2 


(30) 


— occur representing two possible oscillatory 
x ! 


A typical graphical diagram of the curve N[— ,<jd| for one value 

\ X 1 / 

of (i) as varies from 0 to 1 is shown in figure 4. As shown, 

A 1 

for a common frequency the curve for - i intersects the curve 

for G(jco) at two points b and c as demonstrated by equation (30). 


It should be noted in figure 4 that the curves for descending and 
ascending gain have been separated along their abscissa for clarity. 
Actually the curves, because of constant phase angle, are coincident 
and are hereinafter plotted in that manner. 


As shown in reference 1, it is necessary to determine the incre- 
mental increase or decrease in the amplitude of the output Y for an 
incremental increase in X]_ to determine whether the oscillatory con- 
dition corresponds to a stable or unstable limit cycle. Applying this 
analysis shows the first intersection for increasing X-^ to be unstable 

and the second intersection to be stable. Thus, the interpretation of 

figure 4 is that for the region from a to b a stable system 

X 1 

exists in which the amplitude will die out, the region from b to c 
corresponds to an unstable region in which an oscillation will build 
up to the amplitude X 1 at c, and the region from c to d cor- 
responds to a stable region in which the amplitude will reduce to the 
amplitude X-^ at c. Thus, the solutions of equation (30) correspond 
to the stable and unstable points, respectively. 

If in equation ( 30 ) 


> — 

4 


(31) 


a limit-cycle condition will not exist at this frequency. This rela- 
tionship serves to predict the condition for which the net gain will 



always be less than unity. Such a condition results in the variation 
of gain with phase angle for the case in which the curves for the linear 
and nonlinear expressions have no intersection at a common frequency. 

If this is the case, it is desirable to determine the gain margin (the 
difference between the gains of the linear and nonlinear expressions) 
which would be necessary to cause a limit-cycle condition. Therefore, 


the value of the input amplitude for minimum — (or maximum nonlinear 

N -> 


gain) must be determined. 


This is accomplished by differentiating — 


with respect to X-j_ and equating the resulting expression to zero, or 



Rearranging equation (32) and applying trigonometric identities. 



Solving equation (33) for 


+ 2d(x 1 2 - 

D 

— • gives 
A 1 


2\ X /2 

D ) tan T] = 0 


(33) 


A = (L± 
X 1 " \ 


sin 


4 


1/2 


( 3*0 


or the value of for minimum 1 

X-j_ N 

with equal turn-on and turn-off time, 
equals 0 . 71 * 


Thus, for example, for a valve 

that is, t| equals zero, -5- 

X 1 


Another special condition exists, as illustrated in figure 3(a), 
when the valve lag time is longer than the time for the input to turn 
on and off the electronic switch. For this condition when 


(35) - 


2 Oq ^ CDT 

or when substituting from equation ( 19 ) 
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no output appears and the gain 



goes to zero. 


This represents 


a condition on the curve where no physical operation is possible. Thus, 

the maximum possible can be calculated as a function of the valve 

X 1 


switching time. 


The duty cycle , the time percentage of one cycle of operation that 
the reaction motors are on ; can also be determined by the describing 
function. From figure 3(b) , 


Duty cycle 


201 


Substituting equation (20) into equation ( 37 ) gives 


Duty cycle 


2 


cos 


-1 


+ (JDT (k 

A 1 


It 


1) 


( 37 ) 


( 38 ) 


with — and cd determined from the stable equilibrium locus of the 
X 1 

the curves G(jtt)) and - — . 


Case II. Describing Function for Reaction Control and 

Main Thrust Misalinement (Missile Thrusting) 

In this analysis the gain of the secondary system composed of the 
nonlinear components, is determined by imposing on the system an input 
of the form 

X = Xq + cos or r (39) 

In this case, however, the operation may not be symmetrical about the 
abscissa so that the output requires not only a fundamental term but 
also a constant term. Thus, 

Yi = Aq + cos ayr 

The system operation is shown in figure 5 with the output now being 
composed of a summation of the control and misalinement moments, or 
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f(t) - Fjhj - F m b m 

(0 ^ cot a i')' > 


f(t) = -F m b m 

^a-^' < cot < 

> (40) 

f(t) = Fjbj - F m b m 

^ cot •g 


Evaluating the angles gives 

= ” CXq + GOT 


(41) 

a^' = ccq + kGOT 


(42) 


The phase shift a*-, is given by 

= ^( k + 1) = 0 (43) 


As in the previous case a new axis is chosen about giving 


*0 - J J^( F j b j - F m b m)a(“*) - f /" F m b m 


Fj b 1 

- -TT 1 H * Vm 


(tt> 


where 


n _ ^ ( _CL l) _ _ , COT 

o a 0 + 2 ^ ~ ^ 


(45) 


and 


D - X 0 

a~ = cos 

0 X, 


( 46 ) 


Substituting equations (45) and (46) into equation (44) gives 
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F.b, r , d - x A n 

A ° = “XT 2 + 2 (k ' 1] J ' Fmbm (47) 

However, for a steady- state condition to exist A o must equal zero 
giving 



F mV 

F J b i 


m 


which is one condition which must be fulfilled in the case of thrust 
misalinement. 

Evaluating the fundamental or A^_ term gives 



2 r % 

F m^m) cos Jp F m^m cos ^ d(^) 


= * F j b i sin p i 


(49) 


Again substituting equations (45) and (46) gives 


An = 


2 r 

— Fjbj sin cos 


- 1 4^1° + m* 




(50) 


Substituting equations (48) into equation (50), 


2F j b j 

At = — sin 

if 


F mV^ 
\ F J b J , 


(51) 


Rearranging equation (5l) and letting 


F mV 

F j b j 


= R 


(52) 
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gives 


■A-l - 


sin 


M 


(55) 


The nonlinear expression is thus 


w _ ^ _ 2F m b m sin [ R ] -j^( k+ D 

X-l " X-l R e 


(54) 


Substituting equation ( 54 ) into equation (5) yields 

— — = G(Joj) (55) 

2F m b m 8ln [ R ] - ;, f <k+l) 

X x R 6 


A graphical solution is necessary to determine the value of eo 
and X x which satisfy the relationship. Again noting that the phase 

angle of nonlinear expression varies only with 00 a value of to can 
be determined for which the phase angles of the "two expressions are 
equal. Substituting the value of linear gain which corresponds to the 
particular value of oo into equation (55) yields 


2Fmbm sin _[H] 
X ± R 


or 



2 w sin [ R ] 

R 


(56) 


Thus, only one value of results for a limit-cycle condition. Equa- 

tion ( 54 ) shows that the nonlinear gain varies inversely with Xp and 

disturbances in the input amplitude would cause oscillations either to 
build up or to decay back to the value given by equation (5 6). Thus, 
equation (56) represents a stable equilibrium condition. 
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The steady- state error Xq can also "be determined for the equilib- 
rium condition. By equation (48) 


cos 


■1 ^1^0 + £L(k _ i) = 


^1 




or 


Xq = -X^ COS 


fd- k) 

t F j b j 2 


+ D 


(57) 


The duty cycle can be expressed by 


Duty cycle = 


(58) 


or, by equation (4-5) and equation (46), 

„ D - X 


COS 


0 . 0)T 


+ ^(k - 1) 


Duty cycle = 


(59) 


Substituting equation (48) yields 


Duty cycle = 


F m b m 

F j b j 


(60) 


Thus the duty cycle is simply the ratio of the misalinement to the 
reaction-control moment. 

This study does not consider another mode of operation which might 
exist for a condition of smaller main-thrust misalinement. In looking 
again at figure 5 , if the misalinement thrust is sufficiently small com- 
pared with the control thrust, the error signal will cross the dead zone 
-D in the negative half cycle and cause a control- thrust pulse. Such 
a condition is amenable to the describing- function analysis in the same 
manner as the two cases previously discussed. However, it was not con- 
sidered to be as critical for design purposes as the cases investigated. 
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RESULTS AND DISCUSSION 


The results obtained by applying the describing- function technique 
are compared with the results obtained in a simulation of an actual 
vehicle reaction- control system. In the simulation study the actual 
hardware, including reaction motors, was used in conjunction with an 
analog computer which simulated the dynamics of the vehicle. 

The system parameters used in the analytical study are the best 
estimates of the hardware characteristics from the simulation study in 
order that a comparison of the results of the two studies might be made. 
These parameters are presented in table I. In the simulation study, 
exponential functions were used to approximate the peak, nominal, and 
decaying thrust of the main engine. For the analytical study, however, 
only the nominal thrust value was used with an engine burning time deter- 
mined by the value of total impulse used in the simulation study. 


Missile Coasting Phase 

By using the analytical techniques developed in the previous sec- 
tions and the hardware parameters shown in table I, an analysis of the 
coasting phase of flight has been made. Figures 6 and 7 show plots of 
equation (23) for lag times of "slow" and "fast" valves, respectively, 
and a nominal thrust level of 486 pounds. 

Figure 6 shows gain plotted against phase angle for various values 
of frequency and input amplitude for the nonlinear expression 

As noted previously, a single curve can be used to describe G(jcn) with 
frequency as a parameter. However, separate curves are required for 

each value of co for the nonlinear expression - i(jo>). Each of these 

curves is vertical and extends from a gain of infinity to a minimum at 

approximately 27db and then back to a maximum as — varies from 0 to 1. 

X 1 

Since the range of general interest will be between of 0. 1 to 0. 

X 1 

only this range has been shown on most of the curves. Attention should 
again be called to figure 4 to note that this line is double valued. 

As can be seen in figures 6 and 7> increasing frequencies cause the phase 

of to vary from -178° at o> = 0. 25 to -69 0 at 0) = 15* 

X 1 

In figure 6 an intersection of the curves for G(joo) and - ~ 

N 

occurs for a common frequency of 8.14 radians/sec. This corresponds 
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to the frequency where limit-cycle oscillations can take place. At this 

point rp- values of 0. 5^5 and 0..775 occur. As previously discussed with 
X 1 

figure 4, there are two possible modes of stable operation; the input 
amplitude will die out if the system is operating at a rp- value greater 


than 0. 775, and a limit cycle at the amplitude 


■ = D 
1 0. 565 


will occur 


for a — value less than 0.775* Thus a shock-excited limit cycle is 
X 1 

possible, and indeed in the hardware test one did exist at a measured 
frequency of 6.8 radians/sec. The duty cycle predicted from the analysis 
using equation (58) is 58 percent which compares favorably with the 
48-percent duty cycle obtained from the actual tests. 


It could be reasoned that with proper adjustment a — value 

X 1 


greater than 0. 775 could be maintained, and thus no limit cycle would 
exist. Theoretically this may be possible, but in the analysis it was 

shown that a maximum value of also exists. This has also been 

X 1 

plotted as the dotted line on figure 6. Thus, the only permissible 
range, as illustrated by equation (30), in which a l imi t- cycle condi- 


tion will not exist is from 


_D_ 

X 1 


of 0.844 to 0.775* 


It can be argued 


intuitively that this is not a sufficient gain range (approximately 
ldb) for operation without the existence of a limit cycle and that the 
impulse thrust from a reaction motor firing will have sufficient fre- 
quency content at the resonant frequency to cause the system to be shock 
excited into a limit- cycle oscillation. 


Figure 7 shows the same system operation with fast valves. Here 
no intersection of the G(jcn) and - i curves takes place for a com- 
mon frequency under the assumed conditions. Therefore, no limit-cycle 
frequency will exist at this point. This was found to be true in the 
hardware tests. It should be noted that this system has a 6db gain mar- 
gin before limit-cycle oscillations at 17. 1 radians/sec would take place. 

Only one value of rate gain, = 0. 4, was used in the hardware 

tests. A survey of the feedback stabilization gains was made analyti- 
cally in order to determine the best possible gain, K^, for the slow 

values and a nominal thrust level of 486 pounds. The results are shown 
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in figure 8 where the grid of values plotted shows the G( jcu) loci 
for frequencies from 3 to 10 radians/sec and values of Kq from 0.1 

to 1.0 deg/deg/sec. Plotted also are the - i curves for the same 

N 

range of frequencies and the maximum value of the nonlinear gain, N may , 
^corresponding to minimum obtained by using equation (3M» The 
dashed . curve shows the loci of the intersections of the curves G(ju>) 

and - 7 or the limit-cycle conditions. Under the conditions assumed 

N 

the best gain range is for 0. 22 < Kq < 0 . 37 since for these conditions 
the curves do not intersect since the points are below line. 

However, since the gain margin is on the order of ldb it is expected 
that the system operation would be critically dependent upon hardware 
variations. 

If the higher thrust level of 586 pounds was used, a limit-cycle 
condition would exist for all values of rate gain. However, it would 
still be desirable to specify those rate gains which give minimum duty 
cycle since they correspond directly to minimum fuel expenditure. This 
has been determined by use of equations ( 30 ) and ( 38 ) and is plotted in 
figure 9 . For these conditions gains of > 0. 28 give the better 

conditions for minimum duty cycle. 


Missile Burning Phase 

An analysis of the operating conditions for the missile during its 
burning phase has been made by using the previously developed techniques 
and the results compared with the actual analog tests. A graphical 
determination of the limit- cycle conditions using equation ( 55 ) for the 
fast and slow valves is shown in figures 10 and 11 , respectively, at a 
point in time 16 seconds after the ignition of the main engine. 

These figures show that an oscillatory condition will always exist 
at the frequency of intersection. By using equation (60) and the values 
from figure 10 , the duty-cycle and limit-cycle conditions plotted against 
thrust misalinement are shown for the fast valves in figure 12. Also 
plotted in figure 12 are the data from the analog tests. A similar set 
of curves are shown In figure 13 for the operation of the slow valves. 

The figures show a reasonable and conservative estimate for the 
. duty cycle. The calculated values of duty cycle are greater than the - 
corresponding test values. This discrepancy can be accounted for, how- 
ever, by the fact that nominal thrust levels were assumed for the cal- 
culated values. During high-frequency operation of the rocket motors, 
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however, the transient thrust period becomes predominant, and the average 
thrust over the time interval of operation is greater than the nominal 
value. The two results will give identical values in the determination 
of fuel consumption during a specified time interval if average thrust 
values are used in the analytical equation. 

However, for the limit-cycle condition, agreement is not as good 
because the analytical results predict that the frequency of limit cycle 
does not change with thrust misalinement while the test data show a 
definite increase with increasing misalinement. One factor which has 
been neglected in the approximation and which might account for this 
discrepancy is the effect of the higher harmonics on the input to the 
system. Since the amount of harmonic content present is altered for 
different values of R, this could, in turn, alter the limit-cycle fre- 
quency. ~The method for analyzing the effect of higher harmonics is shown 
in reference 2 , but since this method is long and tedious, no attempt has 
been made to apply it to these results. It is true, however, that the 
results do afford a good prediction to the average frequency of the limit 
cycle, as is shown in figure 14. 

The describing function by using equation (57) can also be used to 
predict the steady- state error about which the system oscillates. Fig- 
ure 15 shows that very good agreement exists between the analytical and 
test results. 


Reaction-Control Fuel Consumption 

The values of duty cycle previously determined were used to compute 
the reaction-control fuel consumption for a 50 . 85 - second period of main 
rocket motor burning followed by a 5~ second coast period. The following 
assumptions were used in making the calculations: 

(1) The reaction motor burns for 2 seconds during the first 
5 seconds to overcome initial transients. 

(2) A limit-cycle condition exists for the remainder of the main 
motor burning period. 

( 5 ) A limit-cycle condition exists for the entire 5 ~ second coast 
period when the slow valves are employed. 

(4) When fast valves are used, an arbitrary duty cycle of 10 per- 
cent is assumed during the coast period. 

By using these assumptions, the fuel consumed over the time interval 
was determined by dividing the total impulse by the specific impulse 
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of the fuel. The analytical results are compared with the fuel- consumption 
test results obtained by determining the actual fuel weight differential 
for the test run. The comparison is shown in figure l6. As indicated by 
the figure, excellent agreement exists between the results of the two 
studies. 


CONCLUDING REMARKS 


The analysis, discussion, and comparison with actual hardware tests 
have illustrated that by use of the describing- function technique, it is 
possible to predict by analytical techniques some of the important fac- 
tors involved with the design of an on-off reaction- control system. The 
following conclusions can be stated: 

1. The technique is applicable to the condition of main engine 
thrusting and coasting phases of flight. 

2. The limit-cycle frequency can be estimated for both phases. 

3. The duty cycles can be obtained for both phases. From the duty 
cycle it is possible to estimate the quantity of fuel used. 

k. Hie parameters which affect the existence of a limit-cycle fre- 
quency for the coast phase can be studied. The effect of feedback, 
thrust level, dead band, and valve lag can all be taken into account in 
the design analysis. 

5. The results of the analytical study show good correlation with 
results obtained in an analog simulation including reaction-control 
hardware . 


Langley Research Center, 

National Aeronautics and Space A dmi nistration, 
Langley Field, Va., October 31, i 960 . 
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TABLE I.- SYSTEM PARAMETERS 


Reaction motor characteristics (dynamic response of eq. (2)): 


Nominal thrust level, (Fj) sg > lb 486 and 586 

Damping, £ 0.6 

Undamped natural frequency, o^, radians/sec 45 


Valve response: 

Slow valves - 

Amount of time to turn on, sec . 

Amount of time to turn off, sec 

Fast valves - 

Amount of time to turn on, sec 

Amount of time to turn off, sec 

System dead zone, radian . . 

Main-thrust-misalinement angle, deg 

Main- engine thrust, lb 

Burning time, sec 

Vehicle coast time, sec 

Pitching moment of inertia (0 % t ^ 30), I, 

slug- ft 2 

Reaction- control- thrust moment arm 

(0 £ t $ 30), bj, ft 

Misalinement- thrust moment arm (0 ^ t % 30), 

bm> ft 

Position gain. Kg, deg/deg 

Rate gain, K a, 

« sec 


0. 14 

0. 12 

0. 06 

0.05 

±0. 014 

0. 1, 0. 1769* and - o. 25 

63,000 

30.85 

5 

32,700 - 265t - 5.6t 2 

. . . . 14 + 0.0064t 2 

. . . 15.6 + 0.0064t 2 

1 

0. 4 


lb- sec 


Specific impulse of hydrogen peroxide. 


lb 
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(b) Sequential operation of system. 
Figure 1.- Concluded. 







2 6 



Figure 3.- Operation of nonlinear system (missile 




f (t) 
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cal output along time axis. 





Figure 4.- Limit-cycle stability analysis 
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Figure 5-- Operation of nonlinear system (missile thrusting with main- thrust misalinement) . 




Phase angle, deg 


Figure 6.- Variation of gain with phase angle for limit- cycle determina 1 
tion during period of missile coasting. =0.4 fg C 

t = 0. 14 sec: k = 0 . 85 ; Fi = 486 lb; £ =0.6; = 45 radians/sec. 
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Figure 8.- Variation of limit- cycle frequency with rate gain during 

period of missile coasting. Kn = 0. 4 T = 0. 14 sec; 

^ sec 

k = 0 . 85 ; Fj = 486 lb; £ =0.6; cDq = 45 radians/sec. 





Phase angle, 0 , deg 


Figure 10.- Variation of gain with phase angle for limit-cycle determi 
nation during period of main-engine burning (t = l6 sec). 

=0.4 ~"ggg 6 ^ ) t = 0. 14 sec; k = 0.85; Fj = 486 lb; £ =0.6; 

% = 45 radians/sec. 





Phase angle, 0 , deg 


Figure 11.- Variation of gain with phase angle of limit- cycle determi 
nation during period of main- engine burning (t = 16 sec). 

Kg =0.4 ~ e ~ fgg S ~ > T = 0.0 6 sec; k = 0 . 83 ; Fj = 486 lb; £ =0.6; 
o>n = 45 radians/sec. 






Main-thrust misalinement, deg 


Figure 12.- Variation of limit- cycle frequency and duty cycle with main 
thrust misalinement. (t = 0.06 sec.) 
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Figure 14. - Variation of limit-cycle frequency with valve- response time. 
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Figure l6.- Variation of fuel consumption with main- thrust misalinement. 

(30.85 sec burning, 5 sec coast. ) 
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